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Abstract. We study dynamic heterogeneities in the out-of-equilibrium coarsening 
dynamics of the spherical ferromagnet after a quench from infinite temperature to 
its critical point. A standard way of probing such heterogeneities is by monitoring 
the fluctuations of correlation and susceptibility, coarse-grained over mesoscopic 
regions. We discuss how to define such fluctuating coarse-grained correlations and 
susceptibilities in models where no quenched disorder is present. Our focus for 
the spherical model is on coarse- graining over the whole volume of N spins, which 
requires accounting for 0{N^ 1 ^ 2 ) non-Gaussian fluctuations of the spin variables. The 
latter are treated as a perturbation about the leading order Gaussian statistics. We 
obtain exact results for these quantities, which enable us to characterize the joint 
distribution of correlation and susceptibility fluctuations. We find that this distribution 
is qualitatively different, even for equilibrium above criticality, from the spin-glass 
scenario where correlation and susceptibility fluctuations are linked in a manner akin to 
the fluctuation-dissipation relation between the average correlation and susceptibility. 
Our results show that coarsening at criticality is clearly heterogeneous above the upper 
critical dimension and suggest that, as in other glassy systems, there is a well-defined 
timescale on which fluctuations across thermal histories are largest. Surprisingly, 
however, neither this timescale nor the amplitude of the heterogeneities increase with 
the age of the system, as would be expected from the growing correlation length. Below 
the upper critical dimension, the strength of correlation and susceptibility fluctuations 
varies on a timescale proportional to the age of the system; the corresponding amplitude 
also grows with age, but does not scale with the correlation volume as might have been 
expected naively. 



1. Introduction 

Dynamic heterogeneities, where different local regions of a system evolve on different 
timescales, arise in many non-equilibrium situations. Their existence in supercooled 
liquids and other glassy systems has been probed experimentally using techniques 
including light scattering and confocal microscopy [H El El HI El E] and this has been 
complemented by results from simulation and theory [3 El El OH [HI [121 H31 El IT5] . 
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In systems with quenched disorder such as spin glasses, the disorder itself provides 
an obvious source of heterogeneous dynamics: spins compelled to assume particular local 
configurations by the disorder decorrelate slowly, while less constrained ones can lose 
memory of their configuration at some initial time very quickly. In the aging dynamics 
of such systems after a quench to low temperature, it has been argued that an invariance 
of the global dynamics under reparametrization of time also dominates the fluctuations 
of the local dynamics, with different regions of the system effectively having different 
ages [16j H3, [TU HHj- In order to study such behaviour it is natural to consider two-time 
correlation functions and susceptibilities. For a lattice system with N sites labelled by % 
and spins (or more generally local order parameters) Si these functions are, after spatial 
coarse-graining over the entire finite-sized system^, 

£(Mw) = ^J2c u (t,t w ), ciifto = s t (t)s t (t w ) (i.i) 
x(t,t w ) = — Xu{t, t w ), Xii(^)^w) = Qfo.'^j- y ■ (i # 2) 

Here hi is the field conjugate to Si and is assumed to have been switched on at the waiting 
time t w (measured from the time of preparation of the system, e.g. by quenching), with 
the response measured at the later time t. We have used hats to emphasize that the 
correlator and susceptibility defined above fluctuate across thermal histories (including 
variability in initial conditions); we return below to what this implies for measuring 
X- Averaging over the thermal fluctuations gives the conventional correlation and 
susceptibility, C = (C) and x — (%}• I n a g m g systems these are related by an out- 
of-equilibrium fluctuation-dissipation (FD) relation [20] 



- d twX (t, t w ) = R(t, t w ) = X&b*ld tw C{t, t w ) (1.3) 
where R(t, t w ) is the impulse response function, linked to the susceptibility by 

x(Mw) = f dt'R{t,t') (1.4) 

and X(t,t w ) is the fluctuation-dissipation ratio (FDR). The FDR can be read off from 
the negative slope of a parametric FD plot showing Tx versus C, at fixed t and with t w 
varying along the curve. In equilibrium, the FD theorem (see e.g. [21]) implies that the 
FD plot is a straight line with slope —1, corresponding to X = 1. Out of equilibrium, 
the prediction from local time reparametrization invariance |T6], [T7J [HI HH] is that the 
contour lines of the joint probability distribution of the fluctuating quantities C and 
Tx, at fixed t w and t, follow the local slope of the FD plot constructed from the average 
quantities. 

An obvious question to ask is whether the fluctuations of correlation and 
response obey similar constraints in systems without quenched disorder, where dynamic 

§ In C one should, in principle, subtract a term m(t)m(t w ) with m(t) = (1/N) Si(t), in order to 
get a connected correlator. However, this subtracted term is 0(1/ N) in the scenarios we consider and 
so contributes negligibly both to the 0(1) average and the leading C(iV -1 / 2 ) fluctuation of C. 
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heterogeneities are "self-generated". There is some evidence for an affirmative answer 
from simulations of kinetically constrained models of (structural) glasses [22]. Our 
interest here is in simpler models displaying aging, where we can hope to make progress 
by analytical calculation, namely coarsening systems [231 [23]- It nas been argued that 
in these full time-reparametrization invariance no longer holds, with only time rescaling 
remaining as a symmetry in the long-time aging dynamics [25]. This means that there is 
no obvious reason a priori for the presence of any constraint linking local correlation and 
response fluctuations to the average FD relation, and our main aim will be to investigate 
what effects this has on the distributions of the fluctuating quantities. 

Coarsening systems are magnetic systems - or their analogues in gas-liquid phase 
separation, demixing of binary liquids etc - quenched from the high-temperature phase 
to or below their critical temperature, T c , (see e.g. [26] [27] and the review [28]). During 
the phase ordering (below T c ) or the critical relaxation (at T c ), aging occurs due to 
the growth of a length scale (domain size or correlation length respectively) [29], and 
in an infinitely large system equilibrium is never reached. At any time there are 
spins at the interfaces that behave quite differently from the ones inside a domain, 
so that the origin of dynamic heterogeneities and the associated length scale have a 
clear physical interpretation. As before one can then consider the resulting fluctuations 
around average FD plots. For critical coarsening the analysis has to be modified 
slightly: here the FD plots can in fact hide the interesting aging behaviour of the 
FDR J^sT [261 [271 [301 EH [321 [331 [33] . The latter typically is a smooth function of t w /t [26] , 
so that also the fluctuation effects have to be analysed in terms of the same scaling 
variable. 

In this paper we study the fluctuations in the Langevin dynamics of finite-size 
spherical ferromagnets [351 [36J after a quench from equilibrium at infinite temperature to 
some low temperature T. We will focus mainly on quenches to the critical temperature, 
but comment also on the behaviour in the equilibrium region above. Our calculation 
is based on a leading order expansion in 1/N of the non-Gaussian fluctuations of the 
spins, so that we are effectively considering systems of finite but large size N. The 
nature of the expansion prevents us from accessing the phase ordering below T c , where 
as soon as domains are formed the non-Gaussian fluctuations become dominant rather 
than a small perturbation about the leading order statistics, which are Gaussian in the 
spherical model. Our results therefore complement those of Ref. |25j, where for the 
case of zero temperature and N — > oo the fluctuations of correlations were analysed, for 
coarse-graining volumes ranging from a single site to much larger than the correlation 
length. The overall scaling of correlation fluctuations in coarsening below T c is well 
understood and has a clear interpretation in terms of the growing domain size [23]. Our 
study fills an important gap in allowing us to look at coarsening at criticality, where the 
connection between dynamical heterogeneities and a growing correlation length is much 
less clear. 

We will first describe, in Sec. [2} how we characterize the fluctuations of coarse- 
grained correlation and susceptibilities. A discussion is then given of possible alternative 
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definitions of these fluctuating quantities. Some of these can be ruled out as less useful 
because they give different scalings with system size for the variance of correlations and 
susceptibilities. 

Our analysis proper starts in Sec. [3] with the derivation of the correlation and 
susceptibility variances and the covariance. We will give general, exact expressions for 
these quantities in terms of three-time kernels, D and D x , which will constitute the basis 
for all our further analysis. Our interest will be in the out-of-equilibrium dynamics of 
the system after a quench from an initial state of equilibrium at high temperature to its 
critical temperature. Before considering non-equilibrium, though, we will analyse briefly 
the situation of a quench to above criticality, (in Sec. HI), where the equilibration process 
is fast and a genuine equilibrium dynamics takes place. We will first derive the general 
equilibrium expressions of the relevant quantities and later specilize to the case of high 
temperature. Even here, the results are new as far we know. In Sec. |5]we will turn to the 
more interesting case of quenches to criticality. In the regime of small time differences 
the critical dynamics displayed by the system is stationary and one can look at the 
equilibrium situation. For larger time differences the aymptotic dynamics in d > 4 is 
essentially given by the equilibrium one modulated by relatively weak aging corrections, 
whereas for d < 4 one needs to look directly at the non-equilibrium situation. The latter 
is discussed in Sec. |6j We summarize and look forward to avenues for future research in 
Seed 



2. Definition of fluctuating correlation and susceptibility 

In this section we describe first how we will characterize the fluctuations of the correlator 
and susceptibility defined in Eqs. (11.11)1.21) . We then discuss, and largely rule out, 
alternative ways of defining fluctuating correlation and response functions that are 
coarse-grained across an entire finite system. 

To the leading order in 1/N, i.e. inverse system size, that we keep, the joint 
distribution of the fluctuating correlation and susceptibility is Gaussian and therefore 
fully characterized by its second moments. Specifically, denoting the fluctuations 
SC = C — C and 5\ = X ~ X-, we wm study the behaviour of the variances 

V c (t,t w ) =N([5C(t,t w )} 2 ) (2.1) 

V x (t,t w ) = NT 2 ([S X (t,t w )} 2 ) (2.2) 

and of the covariance 

V Cx (t, t w ) = NT (6C(t, t w )5 X (t, t w )) (2.3) 

Factors of N have been included here to make all three quantities of order unity. Also, as 
suggested by the equilibrium FD theorem, we have scaled the susceptibility fluctuation 
Sx by a factor T to obtain a quantity with the same dimension as SC. 

Writing out the variance of the fluctuations of the correlation explicitly as 

V c (t,t w ) = ^E[<^(Mw)<^(Mw)> - (C«(t,0> O)] (2-4) 
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one sees that this is none other than the by now standard four-point correlation 
function used to characterize heterogeneous dynamics, often denoted C 4 or xa [TTj . For 
coarsening below T c in spatial dimension d, the amplitude of this quantity scales with 
~ t d J 2 [23|, where f(t w ) ~ # 2 is the growing domain size. At criticality, £(£ w ) 

- now the lengthscale of regions across which equilibrium correlations are established 

- will still grow with the same exponent but it is less obvious how it enters Vq and 
the corresponding susceptibility fluctuations. Our explicit results will shed light on this 
question. 

There are other fluctuating correlation and response functions we could have 
considered. Firstly, instead of coarse-graining over the entire finite-sized system as 
in (II. II) and (11.21) . we could have coarse-grained over regions of finite size (and then 
gathered statistics also across all possible centre points of such regions). However, 
for the spherical ferromagnet such a locally coarse-grained susceptibility has negligible 
fluctuations compared to those of the correlation, as demonstrated (in the context of the 
leading order Gaussian spin statistics) in [25]. This is why we focus on global coarse- 
graining, for which non-Gaussian effects make also the response fluctuations non-trivial. 

Secondly, our x from ( 11.21) requires that we measure separately, for a given noise 
history, all local susceptibilities Xu- For Langevin dynamics as studied in this paper this 
does not present a problem since the differentiation w.r.t. the local field hi can be carried 
out directly, and an explicit equation of motion for xu be written down, in the spirit 
of a slave estimator (see e.g. [37J). However, already for Markov dynamics simulated 
via a Monte Carlo scheme it becomes necessary in principle to rerun the dynamics A" 
times, each time switching on one of the local fields, unless specially crafted "field-free" 
methods are used [381 HO] . It is tempting to avoid this difficulty by using a standard 
trick for obtaining local responses [U]: one could consider the observable A = J2i e iSi, 
with the 6j quenched zero mean random variables. The response of A to its conjugate 
field, scaled by 1/N to give a result of order unity, is then 



u 

It can be measured with a single rerun of the history (although of course even this 
is likely to be impossible in a real rather than a numerical experiment). The above 
procedure, employed in [161 HE], should reduce to the definition ( 11.21) if the random 
field amplitudes e« are drawn without spatial correlation so that, on averaging over 
their distribution as indicated by the overbar, ejej = <5y. One has to bear in mind, 
however, that the randomness in the 6j may induce additional fluctuations in Xe which 
are not present in x- This does not appear to have been the case for the systems studied 
in [HI HE], presumably due to the presence there of quenched disorder. 

In our case, on the other hand, the variance of Xe would be genuinely larger than 
that of x- This can be seen by considering 



N (( S Xe) 2 ) = T7^ e * e i e ^ ( S Xij$Xkl) = J2 ITP^^i+j'^+k^i+l' (5xw 8 Xk>v) (2.6) 





1 



ijkl 
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where we have used that after thermal averaging the response statistics in a system 
without quenched disorder are translationally invariant. As the statistics of the e are 
defined to be likewise translationally invariant, the normalized sum over i can be replaced 
to leading order by a disorder average, giving 

N (( s Xe) 2 ) = ^oeJWW(Sxoj'Sxk'i') = j7 X) e i e 3 e k e i {$Xij$Xki) (2.7) 

j'k'V ijkl 

The fourth-order disorder average gives 

Qejefce/ = ^j tk^i + U^k CjCi + CjCfc (2.8) 

This is exactly true if the e are taken as Gaussian variables; for e.g. binary variables 
€i = ±1 one gets an extra term —2SijSi^8u but this makes a subleading contribution in 
1/N. Overall one has 

N((5Xe) 2 ) = ± E [{Wxjj) + 2((^) 2 )] (2-9) 
y 

Comparing with (11.21) one sees that this is indeed larger than N((5x) 2 ), by the second 
term in the square brackets of ( 12. 91) . As an aside, we note that ( 12.71) can be written as 
((^Xe) 2 ) = ((^Xe) 2 ), the variance of 5\e is self-averaging with respect to the sampling 
of the field amplitudes e. The same argument can be applied to all other moments of 
Xei so that the entire distribution P(Xe) is self-averaging. The increased variance of 
Xe is therefore present for any given sample of the e. It does not arise, as one might 
alternatively have suspected, by Xe for each sample e having a distribution similar to 
that of x but with a shifted mean that fluctuates with e. 

The difference between Xe and x can be avoided by averaging over a sufficiently 
large number of different configurations of the 6j. (This of course means that an 
appropriate number of reruns of each thermal history are required, defeating to a certain 
extent the object of working with the random field amplitudes e$.) One thus effectively 
"preaverages" over e; allowing for a general covariance e^ej = this gives 

x = ^Y. e aXij (2.10) 

It is this form that we will use in the calculations below, with short-ranged (so that 
J2j e ij — The extreme long-range case = 1 corresponds to spatially uniform, 

non-disordered, fields = 1 and so would be easiest to measure, with only a single rerun 
of the thermal history. The observable A then simplifies to A = J2i Si = Nm so that x 
becomes the magnetization susceptibility Xm- Its (scaled) variance can be written as 

N({Sx m ?) = ^Y,^Xki) (2.11) 

ijkl 

Significant contributions to the sum are expected to arise only when all sites i, j, k 
and I are close to each other spatially (i, j must not be too far apart to give a sizable 
response at all, similarly for k, I, and then these two pairs of sites need to be close to 
each other to have correlated response fluctuations), giving an 0(1) result as for the 
other susceptibilities considered so far. 
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The reason why we will not consider Xm further is that the corresponding correlation 
function C m = (1/N) &ij has a variance that is much larger, by a factor of order N. 
To see this, write C m (t,t w ) = Nm(t)m(t w ) and consider the simplest case of the equal- 
time correlation C m (t,t) = \^/Nm{t)] 2 . The magnetization has fluctuations of order 
1/yfN around zero, so that y/~Nm(t) has zero mean and fluctuations of order unity (or, 
more precisely, of order (, d ^ 2 (t) in d dimensions). The correlation function C m (t,t) = 
[\/Nm(t)] 2 is therefore also of order unity but, crucially, has fluctuations of the same 
order. It follows that N ((5C m ) 2 ) is of order N as claimed. The same argument applies 
to C e defined in analogy with (12. 5ft : one writes C t (t,t w ) = (l/N)J2ij e i e jCij(t,t w ) = 
N[A(t)/N][A(t w )/N] with the staggered magnetization A(t)/N = (l/N)'£ i € i S i (t) 
which scales in the same way as m(t). 

One can phrase the argument for these large correlation fluctuations differently, to 
see more clearly where the difference to the susceptibility fluctuations arises. Taking 
the magnetization correlator, one has by analogy with (12. lip 

N((6C m ) 2 ) = 1 ^(SdjSCu) (2.12) 

ijkl 

When the sites i, j are far apart, Cjj is small and so SCij = Cij — Cij ~ (7^ = Si(t)Sj(t w ). 
But then SC^SCm ~ Si(t)Sj(t w )Sk(t)Si(t w ) can still be substantial as long as i, k are 
close and similarly j, I (or i, I and j, k). There are 0(N 2 ) such terms in the sum (I2.12p . 
giving a scaled variance of C m of O(N) as claimed. The same argument can be applied 
to the variance of C e for uncorrelated ej, which is given by an expression analogous 

to AM]). 

Only by preaveraging over the field amplitudes does one obtain a correlation 
function with fluctuations of the same order as the corresponding susceptibility. By 
analogy with (12.101) . this correlator can be written as 

C* = ^E^A (2-13) 

ij 

In summary, the only sensible definitions of the fluctuating correlation and 
susceptibility that involve coarse-graining across the entire system appear to be (12.131) 
and (12.10j) ; other definitions involving quenched field amplitudes e« without preaveraging 
lead to correlation variances that are larger than those of the susceptibility by a factor 
of O(N). The arguments we have given apply quite generically for systems without 
quenched disorder, [jj In the spherical model the situation is, in fact, somewhat more 
complicated because of the effective long-range interaction between spins arising from 
the spherical constraint. The resulting weak but long-range correlations lead to extra 
contributions to the fluctuations of C but without changing the scaling with iV; for the 

|| They apply also to coarse-graining over finite volumes, as long as we are considering moderate 
timescales where the typical correlation volume remains much smaller than the coarse-graining volume: 
again the alternative definitions that we have considered would give a correlation variance much larger 
than the susceptibility variance, by a factor of the order the ratio of coarse-graining volume to correlation 
volume. 
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susceptibility, these long-range terms provide the only source of fluctuations but again 
the scaling with N is unaffected. 

We will retain the preaveraged field correlations as essentially arbitrary short- 
ranged quantities during the initial part of our analysis, but then simplify in the concrete 
evaluation to the case of coarse-grained local quantities, = dij, effectively returning 
to the definitions fll.l|ll.2l) given in the introduction. Investigation of the more general 
case could be an interesting subject of future work; indeed, only for zero temperature, 
where spins within domains are fully correlated with each other, would one expect to 
obtain correlation fluctuations equivalent to those for the local case. 



3. Setup of calculation 

We analyse the mesoscopic fluctuations in the dynamics of the spherical ferromagnet 
H=W(S l -S J ) 2 (3.1) 

where the sum runs over all nearest neighbour (n.n.) pairs on a <i-dimensional unitary 
(hyper-)cubic lattice. The spins Si are real variables at each of the N lattice sites r i; 
subject to the spherical constraint Yh Sf — N. 

The Langevin equation for this system subject to thermal noise can be written 

as 1 33 1 



dtSi = {z {t) + N- 1 / 2 z 1 (t))S l (3.2) 

where z (t) is the Lagrange multiplier implementing the spherical constraint and N~ l l 2 zi 
is its leading fluctuation of 0(N~ 1 ^ 2 ). The latter is conventionally neglected in the 
Gaussian theory, and this is justified for observables that probe correlations on scales 
small compared to the size of the system. For globally coarse-grained quantities like 
our C and x, on the other hand, one requires the correlations of all the spins of the 
system. The fluctuations of 0(N~ 1 ' 2 ) are then no longer negligible and the Gaussian 
theory becomes invalid [33] . One can also write f l 3 . 2 H in terms of the discrete (lattice) 
Laplacian Q, which takes the values Qu = 2d on the diagonal and Qij = — 1 for n.n. 
sites 

d t S t = • C " (z (t) + N- 1 / 2 z 1 (t))S i (3.3) 

j 

One expects that the fluctuations in the Lagrange multiplier of 0(N~ l l 2 ) induce non- 
Gaussian fluctuations in the spin variables of the same order. To account for this we 
decompose the spin variables as Si = Sj + N~ 1 / 2 r i , where gives the limiting result for 
iV — > oo, which has purely Gaussian statistics, and N~ 1 ^ 2 r i is the leading-order non- 
Gaussian fluctuation correction. Inserting this decomposition into (13. 3D and collecting 
terms of 0(1) and 0(N~ 1 / 2 ) gives 

d t Si = - ^2 + & ~ z a(t)si (3.4) 
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dtn = - 53 QijTj - z (t)ri - Zi(t)s, (3.5) 
j 

In terms of the Fourier components S" q = J2i Si exp(— iq • r^) of the spins the Gaussian 
dynamics (13.41) reads 

d t S tl = -(u tl + z Q (t))S^ + ^ (3.6) 

where u; q = 2X)f=i(l — cosg a ). Its solution with initial condition at time t w is 

S q (t) = i? q (t, t w )S q (t w ) + /* eft' i2q(t, t')f q (*') (3.7) 
given in terms of the two-time Fourier mode response function 



i? q (Mw) = exp (-u(t - t w ) - jf tft'z(i)) 
where the subscript q in uj n has been omitted and 



9\pv) „-o;(t-t w ) 



g{t)=e^ V (2^dt'z Q {t')\. (3.9) 



The two-time correlator in the Gaussain theory reads C q (t, t w ) = (1/iV) (^Sq(t)Sq(U 
and follows from (13.71) by propagating the equal-time correlator C q (t w , t w ) 
(l/N) ^Sq(t w ) S^(t w )\ from initial time t w to final time i 

Cq(Mw) = R q (t,t w )Cq(t w ,t w ) (3.10) 

Once the function g(t) is known, these results capture all of the leading order Gaussian 
dynamics of the spins. Notice that the impulse reponse (13.81) is deterministic: there are 
no response fluctuations within the Gaussian theory. 

To determine the non-Gaussian corrections (13.51) . one needs to have an expression 
for the Lagrange multiplier fluctuations. This can be worked out as [33] 

Zx = - J dt' L(t,t')A{t') (3.11) 

where A(£) is an 0(1) quantity describing the fluctuations of the squared length of the 
Gaussian spin variables Sj 

A(f) = 4s £(*?(*)-!) (3-12) 
ViV l 

and L is the inverse operator of the kernel K 

K(t,t') = 1 ( Sl (t)R im (t,t')s m (t')) = jjR im (t,t')C im (t,t') (3.13) 

defined by 

J dt'K{t, t')L{t', t w ) = 5{t - t w ) (3.14) 

Both K and L are causal, i.e. they vanish for t < £ w . Above and in what follows 
the summation convention for repeated indices is used. In (13. 13ft . Ri m is the inverse 
Fourier transform of (JHHJ), R im = (1/A^) E q e iq (r ^ rm) i? q , with the sum running over 
the N appropriate wavevectors q; for even L their components are integers in the range 
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—L/2 . . . — 1, 0, 1, . . . L/2 — 1 multiplied by an overall factor 2n/L. When considering 
continuous functions of q this sum can be replaced by the integral / (dq) , where we 
abbreviate (dq) = dq/ (2ir) d , and the integral runs over the first Brillouin zone of 
the hypercubic lattice, i.e. q G [— 7r, n} d ; this simplification will apply throughout our 
analysis. In Fourier space the kernel (13. 13ft then reads 

K(t, t') = J (dq) i? q (t, f )C q (*, (3.15) 

The non-Gaussian corrections to the spins are determined by solving the dynamical 
equation (13. 5p . and can be expressed in terms of the Gaussian spins as 

n(t) = ~J dt'dt"R lk (t,t')s k (t')L(t',t")A(t") (3.16) 

As explained in the introduction, the object of our study are the globally coarse 
grained correlation and susceptibility functions, 

C(t,t w ) = ±'£e iS S i (t)S j (t w ) (3.17) 

ij 

1 dSJt) 

x{t - K)= N^'8hM) (3 ' 18) 



For the correlation function we insert the spin decomposition Si = Sj + r^/V N and 
expand to the order 1/y/N of the fluctuations we are interested in: 



C(t, t w ) — — ^2 e ij 

ij 



Si(t)sj(t w ) + — != (n(t)sj(t w ) + Si(t)r;j(t w )) 



(3.19) 



To obtain the corresponding susceptibility we need to expand the spin variables in 
both the magnetic field and N~ l l 2 . More specifically, consider perturbing the system 
by an external field hi = hei that couples linearly to the spins Si. We keep the q fixed 
initially and perform the preaveraging afterwards. The equation of motion in presence 
of the perturbation reads 

dtSi = -fiyS, - (zo(t) + Zl{t) j^^ j S l + h l + 6 (3.20) 

where now a change in the Lagrange multiplier induced by the field perturbation, 
N^^hAzi, is present in addition to the fluctuating component Z\ of 0(N~ 1 ^ 2 ) of 
the unperturbed dynamics. (One can show that there is no 0(h) perturbation in the 
Lagrange multiplier; such a term appears only if the system has a finite magnetization 
and is perturbed by a uniform field [33].) Inserting the corresponding expansion for the 
spin variables 

c ,7 A , rj + hAn 

Si = Si + hAsi H j^j- 2 — (3.21) 

and collecting the O(N ) terms gives to O(h ) the unperturbed equation of motion for 
Si and to 0(h x ) a deterministic equation for the perturbed components 

d t Asi = -ttijAsj - z (t)Asi + 6i (3.22) 
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Integrated in time with the condition Asj(t) = for t < t w this gives 

As i {t)=x i j{t,t vr )e j (3.23) 

where Xij 1S the non-fluctuating Gaussian susceptibility Xij'(^w) — // w dif R%j. 

Gathering the 0(N~ 1 ^ 2 ) terms in (13.201) . on the other hand, gives to 0(h°) 
equation (13. 5p . as expected, and to 0{h}) a new equation 

d t Ar l = -QijArj - z (t)Ar l - Az 1 (t)s l - z^Si (3.24) 

for the Ari, with solution 

An(t) = - f dt' RijitJ^Az^Sjit') + z^Asjit')}. (3.25) 

The fluctuations in the response thus arise from the fluctuations of the Lagrange 
parameter, as anticipated. The Azi term can be worked out by imposing that, due 
to the spherical constraint, N' 1 J2i S 2 {t) = 1 at all times. Using (13. 21 j) . this implies 
that the quantity 

1 E(^ 2 - 1) = ^ E U - 1 + ^^w 11 + 2hAs > s * + 2hAs *7k) ^ 



A 1 „ 2h^ . 2h „ A 

^ + 2 ^ E + Jf E Si * + ^ E *Ar< 



2h 

+ ^E^ (3-27) 

must vanish to the leading order in h, iV _1//2 and hN' 1 ^ 2 ; we have temporarily re- 
instated the summation signs for clarity. To make progress, let us note that the first 
two terms on the r.h.s. of (I3T2~TD cancel to 0{N~ 1 / 2 )] this is in fact how the corrections 
r« are determined |33j. In the third and fifth term we can insert the deterministic 
quantities Asj from (I3.23p . Since the are "driven" by the Sj according to (I3.16p . they 
will only have spatial correlations of finite range. Thus (1/N) J2i r i^ s i i n the fifth term 
is (9(iV -1 / 2 ), making this contribution 0(/iiV _1 ) overall and subleading compared to the 
third term, which is 0(hN~ 1 / 2 ). So we need to impose 

1 / Ar\ 

^Si Asj + -^ =0 (3.28) 



to 0(1), which yields using ( 13.251) 

— / dt' R im {t,t')s l {t)s m {t')Az 1 {t') = -= Si (t)Asi(t) 

f dt' Rinit^SiWAsnWztf) (3.29) 

In the second term on the RHS, Ri m (t,t')As m (t') = Ri m (t,t')xmj(t',t w )ej is a 
deterministic 0(1) quantity which is then summed over sites i multiplied by the short- 
range correlated Si(t). Together with the 1/N prefactor this gives a negligible 0(N~ 1 / 2 ) 
contribution. On the LHS, (l/N)R im (t,t')si(t)s m (t') has fluctuations of 0(N~ 1 / 2 ) 
which can likewise be neglected compared to its 0(1) average; the latter equals K(t,t') 
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from (13.151) . Inverting the resulting convolution dt' K(t,t')Azi(t') using f!3.14p one 
finds as the solution of f!3.29j) 

1 



A Zl (t') 



dt"L(t',t") Sl (t")A Si (t") 



(3.30) 



With this we can now write down the susceptibility for the given set of ej, as defined 
in (12.51) . Noticing that Xij e j is the response of spin i, given by the O(h) terms from (13. 21 j) . 
one gets 



Xe(^5 



1 






■( 


= N € 




1 




= N e 




X 


\I 



N 



An 



N-^ti f dt' Rim(t,t') X 



dt" s m (t')L(t',t")^=s n (t")As n (t") + As m (t')zi(t') 



(3.31) 

The first term is the non-fluctuating Gaussian contribution. The fluctuating remainder 
becomes, once we insert ( 13.231) and preaverage over the e^, 

rt 



Sx(t,t v 



-N- 2 e 



dt'dt" R im (t, t')s m (t')L(t' , t")s n (t") Xnj (t", tw) 



- N-^e^ f dt'R m (t,t')x mj (t',t w ) Zl (t') 

J tw 



(3.32) 



In this expression the first term is 0(1/N): the sum over i, j gives an 0(1) translation 
invariant function €ijR im (t,t')xnj(t' \tw) , an d s m (t')s n (t") can be replaced by its average 
C m n{t',t") to leading order; with the 1/N 2 prefactor and the summation over m,n 
one gets 0(1/N) overall. (The neglected fluctuations of s m (t')s n (t") will give an even 
smaller correction, of 0(N~ 3 / 2 ).) In the second term one argues similarly that the sum 
of (l/N)eijR im Xmj over i,j,m is 0{1). Since zi is scaled to be C(l), it is then this 
term that provides the leading susceptibility fluctuation of 0(N^ 1 ^ 2 ). Inserting (13.111) 
and (I3.12p we can finally write 
1 

2" 

We have dropped the integration limits since these are enforced automatically by 
causality of Ri m , L and Xmj- 

In order to study the fluctuations of globally coarse-grained quantities around their 
mean values, we will consider their variances and covariance, defined in (12.11) . ( 12. 21 ) 
and (12.31) . For the correlation variance one has, by inserting (I3.19P into (12.11) and 
multiplying out, 



--N- 2 e 



it 



J dt'dt" R im (t, t')L(t\ t") Y,(s 2 n (t") - l) Xm3 (t'\ t w ) (3.33) 



Vc(Mw) = j^tijtki {(sj(t)s j (t w )s fc (t)s i (t w )) / 
i 

(r i (t)s i (t w )s fc (t)s z (t w ))' + (s i (t)r i (t w )s fc (t)s i (t w )) / 



+ 



+ (siitfsj^rkitfsiitw))' + (s i (t)s j (t w )s fe (t)r ; (t w ))' 



+ 



1 r 

N 



(ri^Sjit^rk^s^U))' + (r i (t)s j (t w )s fc (t)r z (t w )) / 
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+ (s i (t)r i (t w )r fe (t)^(tw))' + (s i (*)r i (* w )s fc (*)r,(t w )) / ] } (3.34) 

where the prime on the averages indicates that the corresponding disconnected 

contributions arising from (Cij){Cki) are to be subtracted. The susceptibility 
variance f)2.2p reads 

V x (t,t w ) = ^eyew dt'dt''dt' w dCR l ^t')L(t',t'')R kp (t,t'jL(tl,C) 

X £(4(0 " 1) £( S r(0 - l)Xrn#,*wWOw)) (3.35) 
n r I 

while for the cross correlation (12.31) . one has 

T ft , 

Vcx(Mw) = -j^ e ij e M J t dtdt L(t ,t )R kp (t,t)x P i(t ,U) 

X ( + ^ ri(t) ) ( Sj(tw) + ^ ri(tw) ) ^ (t " } " 1} ) (3 ' 36) 

Since all the quantities appearing in the averages can be expressed, via f)3.16p . in terms 
of Gaussian variables Sj, we can use Wick's theorem to perform the averaging. This gives 
a sum over all possible pairings of the Gaussian variables, each contributing a product 
of correlation functions. In the primed averages in Vc, pairings that do not couple the 
index groups [ij] and [kl] need to be discarded, and similarly in Vc x . Fortunately, many 
other pairings can also be dropped because they give subleading terms in 1/N. We omit 
the details as the reasoning is analogous to that in [33], and state the results only for 
the coarse-grained local correlation and susceptibility (ejj = Sij). 
In order to make the expressions more manageable let us define 

CC(t,t w ) = J(dq)Cl(t,t w ) (3.37) 

and the following three-time function (we use the same symbol as the number of 
arguments will make it clear which function is meant; note that CC(t, t, t') = CC(t, t')) 



CC(t, t w , t') = J (dq) C q (t, OC q (tw, (3.38) 
We also introduce 

Dit^t 1 ) = J dt"L(t",t') J (dq)R (l (t 1 ,t")C cl (t 2 ,t") (3.39) 

as well as 

D(t, t w , t') = X - [D(t, t w , t') + D{t w , t, t')] (3.40) 

which is the symmetrized version of (13.391) . Note that D(t,t w ,t') is causal in the sense 
that it vanishes for t' > t. In terms of these functions the correlation variance takes the 
compact form 

V c (t, K) = J (dq) C q (t, £)C q (t w , t w ) + CC(t, t w ) - 4 J dt'D(t, t w , t')CC(t, £ W) t') 

+ 2 J dt'dt' w D(t, t w , t')D(t, t w , t' w )CC(t', t'J (3.41) 
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Similarly one can define 

£>*(*, t w , t') = T J alt" Lit!', | (d<?) i? q (t, t")x q (t", t w ) (3.42) 
and express the susceptibility variance as 

V x (t, K) = \J dt'dt' w D*{t, t w , t')D*(t, t w , t' w )CC(t', t'J (3.43) 

In the susceptibility the times are already ordered and we do not need to consider a 
symmetrized version of D x ; D x is causal in the same sense as D. The covariance, finally, 
can be expressed in terms of the same functions as 

V Cx (t,t w ) = - J dt'D^t^,t')CC{t,t w ,t') 

+ J dt'dt' w D(t, t w , t')D x (t, t w , t'JCC(t', O (3.44) 

All the properties of the (co)variances can now be obtained from the behaviour of the 
functions CC, D and D x . To understand the general structure of D and D x we first 
recall [33] that L, the inverse kernel of K, has the from 

L(t, t') = 8'(t - t') + 2T5(t - t') - L {2) (t, t') (3.45) 

where L^ 2 \t,t') vanishes for t' > t, has a jump discontinuity at t' = t and is expected 
to be smooth and positive for t' < t. The singular terms are consequences of the fact 
that K(t,t') vanishes for t' > t and has equal-time value and slope 

lim K(t, t') = 1, dt>K(t, t') \ t=t ,+ = 2T (3.46) 

The remaining ingredient in D is the function E(t 1 ,t 2 ,t") = J(dq) Rq(t%, t")Cq(t2, t"). 
This vanishes for t" > t\ because of the causality of i? q , and has a jump of size 
/ (dq) Cq(t2, t\) = C(t 2 , ti) as t" decreases past t±. If t 2 < t\, E actually remains constant 
at this value down to t" = t 2 because the ^''-dependence in Rq(ti,t")Cq(t", t 2 ) cancels as 
a consequence of (13.81) and (13.101) . For t" < min(ti,t 2 ), one can use the same identities 
to express E in terms of the kernel K: the q-dependence (via uj) of Rq(ti, t")Rq(t2, t") is 
the same as that of R\{(t\ + t 2 )/2, t"), and accounting for the remaining proportionality 
factors results in Efa, t 2 ,t") = g((t 1 +t 2 )/2)g- 1/2 (t 1 )g- 1/2 (t 2 )K((t 1 +t 2 )/2,t"). Carrying 
out the ^''-integral in (13.391) and exploiting the decomposition (I3.45P of L then gives for 
D the general form 

D(t, t w , = -C(t, t w ) [6(t - f) + 5(t w - 0] + Mt, tw, t'W ~ ^) 

+ D 2 {t,t w ,t')6{t w -t') (3.47) 
where the continuous pieces for t' above and below t w respectively are 

£i(Mw,*') = \c(t,t w ) (2T- £dt"LW(t",t'^ (3.48) 

and, abbreviating t — (t + t w )/2, 

D 2 (t,t w ,t')= - \c{t,t w ) f dt" m (f, if) + 9{f) / ( (~ + 2t) K(t, t') 
2 Jt - Jg(t)g(t w ){\ dt' J 
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dt"K{iX)L {2) {t\t')\ (3.49) 



= --C(t,t w ) f dt"L^{t",t')+ . 9 ® -. [* dt"K(t,t")LW(t",t')(3.50) 

2 Jtr, \Jg{t)g{U) Jt - 

The last simplification for D2 arises because, from ( 13.141) and ( 13.451) . the terms in curly 
brackets in (13.491) would cancel exactly if the upper integration limit was i. 

For the corresponding function D x for the susceptibility, the g-integral in (13.421) 
can also be simplified by exploiting the link (11.41) between x q and i? q : 

J (dq) R^t, t")x q (t", t w ) = X (t, t w ) - X (t, t") (3.51) 

This holds for t w < t" < t; otherwise the function on the LHS vanishes due to causality. 
Inserting into (I3.42[) and using again (13.451) gives 

D*(t, t w , if) = T X (t, t w )5(t - if) + D\ (t, t w , tf)6{tf - t w ) + D$(t, tw, t')^(tw - if) (3.52) 
with 

T- X D\ (t, t w , = - R(t, f) + 2T[ X (t, t w ) - X (t, t')] 

- f dt" (t", [ X (t, tw) - X(t, t")} (3.53) 

and 

T^D^t^t') = - f dt"L^(t",t')[ X (t,t w ) - X (t,t")} (3.54) 

Note that the expressions above are general and valid for arbitrary quenches, since we 
have not imposed any restrictions on the form of response, correlation or the kernel 
L. They will therefore form the basis for all further analysis of the correlation and 
susceptibility variances Vq and V x and their covariance Vc x - 

In addition to the variances and covariance themselves we will also consider the 
correlation coefficient 

7 = (3-55) 

which lies in the range —1 ... 1; the extreme values correspond to susceptibility and 
correlation fluctuations being fully correlated, i.e. identical up to a scale factor. The 
joint probability distribution of (C, x) can be more fully characterized by its contour 
lines. Due to the Gaussian nature of the distribution (in our leading order approximation 
in 1/N) the contours are ellipes given by 

vTYiZ)-"** (3 - 56) 

These are centred on (SC,S X ) = (0,0), i.e. on the mean values (C,x)- Geometrically, 
it is then natural to define the direction of the dominant fluctuations as the principal 
axis of the ellipse. We define the negative slope of this as in analogy with the 
FDR X which gives the negative slope of the FD plot relating the mean values Tx 



Dynamic heterogeneities in critical coarsening 



16 



and C. If the predictions for spin glasses summarized in the introduction also apply to 
coarsening systems, one would expect X fl to be close to X. Explicitly, by diagonalizing 
the covariance matrix in (13.561) and finding its largest eigenvector one has 

/ 




\ 




+ 4 7 2 



(3.57) 



In accordance with the definition of Xr as the negative slope of the principal axis, it 
always has the opposite sign of the correlation coefficient 7. We note that the definition 
of Xq, unlike that of 7, depends in principle on the relative scaling of the axes of the 
FD plot. The factors of T included in (12.21) and ( 12. 31) correspond to measuring the 
fluctuation slope from contours in the (C, Tx) plane where the equilibrium FDT is a 
line of slope —1. While not unique, this is certainly the most natural choice. 



4. Quenches to T > T c 

In this section we study quenches to above criticality so that, as discussed above, 
equilibrium is considered. In equilibrium the average correlation and susceptibility 
functions are time translation invariant (TTI) and related by the fluctuation-dissipation 
theorem. We can ask whether FDT-like relations also hold for the fluctuations of 
correlation and susceptibility around their typical values, i.e. whether X& is close to 
unity. 



4-1. Equilibrium expressions for D and D x 

In equilibrium, all functions depend only on time differences, so we will write K(t, t w ) = 
K(At), L(t,t w ) = L(At) and so on, with At = t — t w . For the three-time functions we 
will keep the three separate arguments; for CC(t,t w ,t') this helps to avoid confusion 
with the two-time function CC(t,t w ) = CC(At). 

In order to work out the equilibrium expressions for D and D x , we need first the 
various covariance and response functions, as well as the kernel K and its inverse L. The 
Lagrange multiplier z approaches a constant value z cq at equilibrium, corresponding to 
exponential growth g(t) oc exp(2z cq t) of the function (13.91) . One can then show that 
Cq(t,t) = T/(z cq + uj); since the spherical constraint imposes J(dq) C q (t, t) = 1 at all 
times, be found from the condition 

(dq)—^— = l. (4.1) 

z eq + u 

For the moment we will leave the Lagrange multiplier unrestricted, so that the following 
results will be valid for equilibrium at any temperature > T c . (For T < T c one would 
need to account separately for the q = mode which acquires a nonzero expectation 
value proportional to the equilibrium magnetization.) Later we will consider first 
high temperatures, then generic temperatures above criticality, and finally, in the next 
section, T = T c where z eq vanishes. 
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The exponential behaviour of g(t) in equilibrium reduces the Fourier mode response 
and correlation functions (I3.8P and (13.101) to the simple forms 

i2q(At) = e~ iu+Zcci)At (4.2) 
C q (At) = T e -(*+*c q )A* ^ 

These determine the equilibrium form of the kernel K (13.151) as 

K(At) = ! (dq) T e ~ 2{uj+z ^ At (4.4) 

J LU + Z eq 

and the (average) local correlation and response can be expressed in terms of this as 

C(At) = f (dq) e -(c+ %q )Ai = ^Afj ^ 

J LU -\- Z eq 

TR(At) =T J (dq) e -(^+^) A * = ~K'{^ ) (4.6) 

Tx(At) = I (dq) T (l - e- (w+ ^ )At ) = 1 - ) (4.7) 
while the two and three time versions (I3.37P and f!3.38|) of CC become 
CC(At) = [ (dq) ( T ) e - 2 (-+^) Ai (4.8) 

J \LU + Z eq ) 

CC(t,t w ,t>)=9(t'-t w ) f(dq) (—^-) e~^ +z ^ At 

+ 6(t w - t') [ (dq) ( T ) e -2(^e q )(f-*') Ug\ 

J \LU + z eq ) 

= CC(At/2)9(t' - t w ) + CC(t - t')9(t w - f) (4.10) 

Notice that CC(t,t w ,t') is independent of t! in the regime t w < t'(< t). 

Finally we need the inverse kernel L. Combining (13. 14ft and ( 13.451) . we can express 
its Laplace transform L( 2 \s) as 

L^(s) = s + 2T--J— (4.11) 

Keq(s) 

We will require occasionally the integral of over all times, which follows as 

because K(0) = dt K(t) diverges at T = T c for d < 4 (see Eq. ( 15.2ft below). 

Putting everything together, we get from (I3.47H . ( 13.481) and ( 13.501) the explicit 
equilibrium form for D: 

D(t, t w , t') = \K(f) [5(t - t') + 5(t w - 0] + D x (t, t w , t')9(t' - t w ) 

+ D 2 (t,t w ,t')6(t w -t') (4.13) 
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where 

D 1 (t,t w ,t') = \k(%) {lT - f'* dr L®{t) 



and 



D 2 (t,t w ,t') = -~K(f) f \drL^(r) + [' ' dr K(t - t' - r) L^(t) 



t w —t' 



Similarly one has from ( 13.521) . ( 13.531) and ( 13.541) for D y 



D x (t,t w ,t') 
with 



1-K 



At 



- / dt"L (2 \t" -if) 



and 



D*(t,t w ,t') = - / dt"L^(t"-t') 



(4.14) 



(4.15) 



5(t - +D?(t, t w , t')9(t' - t w )+D$(t, t w , t')6(t w - If) (4.16) 



(4.17) 
(4.18) 



4.2. HighT 

Having derived the general equilibrium expression for the functions D and D x 
at arbitrary temperature, we next study their time dependence in the regime of 
temperatures above criticality, T > T c . First we consider briefly the limit of high 
temperatures, where explicit expressions can be obtained. 

For T — » oo, one sees from (14.11) that the Lagrange multiplier needs to scale 
as z cq = T + 0(1) because the frequencies ui are of order unity and independent 
of T. This suggests a series expansion as z eq = T + a + b/T + 0(1/T 2 ), and by 
substituting into ( 14.11) and using J(dq)uj = 2d and J(dq)u 2 = (2d) 2 + 2d one finds 
—a = b = 2d. The time dependence in the equilibrium functions (14.41) . (14.81) and (14.101) 
through the combination (z eq +uj)At is then equal to TAt to leading order. We therefore 
rescale the time difference with temperature as r = TAt and expand all exponentials 
exp[-(z eq + u)At] = exp[-(l + (u - 2d)/T + 2d/T 2 )r] in 1/T. One finds in this 
way K(t) = e~ 2r (l + Adr 2 /T 2 ), to C(l/T 2 ). The C(l/T 2 ) term in K(r) is needed 
to determine the Laplace transform of from (14. lip , because the leading order 
cancels. Inserting K into (14.111) shows that the first non vanishing term in L^ 2 \s) 
is 0(1/T), which transformed back to rescaled time variables yields the 0(1) result 
L^ 2 \r) = 8dexp(— 2r). To use these results in a systematic high-T expansion up to 
0(1/T 2 ) of the correlation and susceptibility (co-)variances, we need to know to which 
order in 1/T the functions that appear need to be expanded. For large T it is convenient 
to rescale D and D x by a factor T in order to work with quantities of order unity. The 

We remark that although 0(1/T 2 ) terms are needed to determine L^ 2 \ the latter has a value of 
0(1) and not 0(1/T 2 ) as we had mistakenly stated in 33J . Fortunately, this error had no effect on the 
calculations in [33 , since the large-T behaviour of was never used explicitly. 
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10 



Figure 1. Correlation and response variance versus the rescaled time r for T = 15. 
Both show a power law increase for small r and an exponential approach to their limit 
value for large r. 



compensating factor XjT is absorbed into the rescaling of the time integrations that 
lead from D and D x to the (co-)variances. One can check that the terms proportional 
to in the rescaled D and D x are smaller than the others by a factor 1/T 2 , because 
they are always obtained by integrating over time. Therefore L^ 2 \t) is only needed to 
(9(1). Expanding all other functions to order 0(1/T 2 ) and inserting into (I3.4ip . (13.431) 
and (13.441) . all integrals can be done explicitly. One obtains, after some lengthy but 
straightforward algebra, 

1 



V c (r) = l-(l + 2T)e 



-2t 



J^2 



2cfe- 4T -4cfe- 2r (l 



2r 2 + r 3 ) 



+ 2d 



(4.19) 



V x (t) = 1 + (3 + 2r)e 



-2r 



4e" T + 



J^2 



+ de- 2r (13 + lOr + 6r 2 + 4r 3 ) 



2 4 
-de~ 4T - -de~ 3T 

3 3 

-cfe- T (ll + 3r 2 ) + -d 
3 3 



(4.20) 



V, 



Cx 



T 



2d 
3T 2 



[2e" 4r - 3e~ 3r + 6(1 



r 2 )e- 2r 



+ (6r - 5)e _ 



(4.21) 



Plots of ([USD , (PB and (Effl) are shown in figures □ and [2] (left) for T = 15. One can 
study the high temperature limit of the correlation and susceptibility variances directly 
by setting the £>(1/T 2 ) corrections in fICTD and ffl~20l to zero. This shows that for high 
T the correlation and susceptibility variances are monotonically increasing funtions of 
r, starting from zero at r = (i.e. At = 0). This is as expected since C(t w ,t w ) = 1 
cannot fluctuate due to the spherical constraint, while x(t w ,t w ) vanishes trivially. An 
expansion for small r shows that the correlation and susceptibility variances increase 
initially as, respectively, Vo(t) = 2r 2 and V x (t) = (2/3)r 3 . These scalings, including 
the prefactors, will also be found at finite temperature (see below). They show that 
there are significant correlations in the time evolution of C and x; if the fluctuations 



Dynamic heterogeneities in critical coarsening 



20 



had independent increments at different times this would lead to a random walk for the 
fluctuations and hence a much more rapid increase of the variances, V ~ r. 

We cannot infer from the 1/T expansion whether the monotonic behaviour in r 
of the variances holds also for finite temperature. However, the fact that already the 
leading 1/T 2 corrections are non- monotonic in r suggests that the overall variation at 
finite T may also be non- monotonic. Indeed, for the correlation variance Vc we will 
see in Sec. 15.11 using different arguments that a non-monotonic dependence on At (or 
equivalently r) occurs at least in d > 4 and for T not too far above T c . 

In the limit T — > oo, both variances approach the constant value 1 exponentially 
fast in t. For the correlation this can be explained relatively simply: as the spins 
Si(t w ) and Si(t) decorrelate at long times and are also uncorrelated in space for large T, 
C = (l/N)J2iSi(t)Si(t w ) becomes a zero mean Gaussian random variable of variance 
1/N. Consistent with this intuition, the dominant contribution to Vc for large At comes 
from the Gaussian fluctuations which are described by the first two terms in (j3.4ip ; in 
fact, only the first term survives for At — > oo. It should be emphasized, however, that 
the high-T limit does not amount to neglecting all non-Gaussian effects. Indeed, the 
Gaussian terms from (13.411) would give the quite incorrect result Vc = 2 for At = 0. 

We next look at the covariance of correlation and susceptibility, and the 
consequences for the correlation coefficient 7 and the fluctuation FDR Eq. (14.211) 
shows that the covariance is 0(1/T 2 ) for any finite r, and it vanishes in the limits of 
both small and large r as, respectively, 

V Cx (r) » -4dr 4 /(3T 2 ) (4.22) 

and Vc x ~ —Adre^ T /T 2 . Plotting the full expression (14.211) (see Fig. [2] left) shows that 
Vc x is negative not just in these two limits but in fact for all r. 

From the above results one can determine the correlation coefficient, as defined 
in (I3.55p . for high T. For r — ► and in the limit of high temperature, one obtains 
directly from the small-r scaling of the (co-)variances that the correlation coefficient goes 
to zero as 7 ~ — t 3 / 2 /T 2 . For the opposite limit r — > 00 of long times, Vc = Vc x = 1 to 
leading order, as explained. This yields 7 ~ Vc x ~ — 4rfre" r /T 2 . A plot of 7 (see Fig. [2] 
right) shows that like Vc x it is negative for all r, and its modulus is smaller than unity 
as it should be. The scaling with 1/T 2 shows that, for high temperatures, correlation 
and susceptibility fluctuations become increasingly less correlated with each other. 

Studying the fluctuation FDR (I3.57P to characterize the joint distribution of 
correlation and susceptibility requires some care. If we first proceed as above, keeping 
r of order unity fixed and taking T — ► 00, then 7 scales with 1/T 2 as we saw earlier and 
so becomes small compared to the other terms of order unity in ( 13. 57ft . We can then 
expand in 7 to get to leading order 

X & = 7 , (4.23) 



f Vc/V x - JVJVc 
For small r, where 7 ~ — r 3 / 2 /T 2 and V x /Vc ~ r, this gives 



r 2 



*fl ~ J* (4-24) 
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Figure 2. Plot of the covariance and the correlation coefficient 7 versus the rescaled 
time r for T = 15. These show a power law increase for small t and an exponential 
decay for large r. Both are negative throughout, and the modulus of the correlation 
coefficient is much less than one, indicating that the fluctuations of correlation and 
response functions are only weakly correlated. 



For large r, the denominator of (14.231) goes to zero even faster than the numerator, and 
we find 

*fl ~ ^ (4-25) 

However, this result must clearly break down when r becomes too large at finite T, 
as Eq. (14.231) was predicated on 7 being small compared to (Vc/V X Y^ 2 — {V x /Vc) 1 ^ 2 ■ 
To understand what happens in this regime, we use that Vc = 1 + 2d/T 2 and 
Vc x = 1 + 7d/(3T 2 ) to leading order for r — > 00 at finite T, where we need to keep 
the 0(1/T 2 ) corrections. Then 7 = —4dre~ T /T 2 in the outer square root of (I3.57P can 
be neglected as smaller than the other term under this root, giving to leading order a 
temperature-independent exponential increase 




(4.26) 



The crossover between the linear and exponential regimes, Eq. (14.251) and Eq. (14.261) 
respectively, can be shown to be due to the competition, for large T and r, between 
0(l/T 2 ) and 0(e~ T ) terms in fl4~19l) and fl420l) . and therefore in <KE7\\ . and takes place 
at t m 2 InT. This is shown in Fig. [3] on the right, along with (on the left) the crossover 
between the quadratic and the linear regimes, Eq. (14.241) and Eq. (14. 25ft . that occurs at 
shorter times. 

With the expression of the (co-)variances to the required orders at hand, we now plot 
the contour lines of the joint distribution of the fluctuating correlation and susceptibility 
(C, T%). The mean values C = (C) and x — (x) can be read off from (14.51) and (14.71) 
as C(r) = exp(— r) = 1 — T%(r) and produce the straight line of slope —1 expected 
from equilibrium FDT. (Note that we do not need to normalize the FD plot because for 
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Figure 3. Left: Plot of the fluctuation FDR Xr versus the rescaled time r, for d = 3 
and T — 15, 20 as shown in the legend. Xa starts off quadratically for small r, then 
becomes linear, its slope increasing with decreasing T in both regimes. Eventually it 
crosses over into a regime of exponential growth, which is not yet visible in the r-range 
shown. Right: Plot of ln(Xfl) versus t, at T — 15, 20 as indicated in the legend. This 
shows the crossover to the regime of T-independcnt exponential increase at large r, 
which takes place at r 2 In T and is represented by the dotted line on the right of 
the graph. 

our local spin correlations the equal time correlator C(t,t) = 1 always.) Fig. H] shows 
contour lines of the fluctuation distributions for a range of different mean values (i.e. 
different r); their centres lie on the straight equilibrium FDT line. For the purposes 
of this graphical illustration, we have aimed to choose a relatively low temperature, 
as otherwise the covariance becomes too small and all fluctuation contours degenerate 
into ellipses oriented along the C and x _axes - We cannot go too low, of course, as 
otherwise truncating the 1/T expansion cannot be justified. The choice T = 10 is a 
reasonable compromise: the 0(1/T 2 ) corrections to the variances are then significantly 
smaller than the leading order terms. However, 0(1/T 3 ) terms - which can be worked 
out along the lines above - are comparable to the (9 (1/T 2 ) contributions, so the results 
shown are not fully quantitative. We consider d = 3 and set the scale of the contours 
by taking N = 50 and unity for the constant on the RHS of (13. 56ft . The relatively small 
value of N was taken only for better visibility; a larger value would simply shrink all 
ellipses uniformly. 

In summary, the fluctuations of correlation and susceptibility are not linked in a 
manner akin to the equilibrium FDT. The dominant fluctuation direction measured 
by Xfl does not lie on the straight line of slope —1 that locates the average quantities. 
Instead, this direction is along the horizontal (correlation) axis for small times and along 
the vertical (susceptibility) axis for large times. 

4-3. Small time differences 

Next we consider the short-time behaviour of the correlation and susceptibility variances 
at some generic temperature T > T c . For the correlation one expects from the high-T 
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c 



Figure 4. Left: Fluctuations of the coarse-grained local correlation C and 
susceptibility T\ in the high temperature limit. The mean values, indicated by the 
centres of the ellipses, lie on the straight equilibrium FDT line. The ellipses themselves 
show contour lines for several different mean values, corresponding to different scaled 
time differences r = TAt; we chose d = 3, T = 10, N = 50 and set the constant on 
the RHS of (|3.56|) to unity. The principal axis of each ellipse is shown, and represents 
the direction of the biggest fluctuations. This lies on the correlation axis for small r 
(Xfl <C 1) but rotates to lie along the susceptibility axis (X& 3> 1) for large r. This 
rotation of the principal axis is continuous but in this graphical representation happens 
very quickly in the top left corner of the plot: this is because it takes place at large 
values of r, i.e. when the mean correlation has already decayed to a very small value. 
The fluctuation contours in this regime are close to circles and would lie essentially on 
top of the leftmost ellipse shown. Right: A zoom of the top left corner of the plot on 
the left. Here only the direction of the principal axis of each ellipse is drawn, to make 
the rotation towards the x-axis clearer. 



result that the leading term for small At will be 0(At 2 ). We therefore expand D as 

[6(t-t') + S(t w -t')} + 6{t'-t w )D 1 + 9(t w -t')D 2 (4.27) 

£>i = T - T 2 At - ^L (2) (0) (t - t') (4.28) 
D 2 = \\2TL^{U - - £ (2) '(tw - t')]At 2 (4.29) 

8 

Here we have used (13.461) ; Di needs to be expanded only to linear order in quantities of 
O(At) (including t — t') because it is integrated over the range t w < t! < t which is itself 
of O(At). Inserting into f!3.4ip and expanding the remaining functions (14. 8 j) and (14. 10)) . 
the various terms involving L^ 2 ' cancel to leading order. In fact, to 0(At 2 ), D 2 only 
contributes to the single integral in (I3.4ip . where the remaining three-time function 
CC is set to its zero-order value, and to the double integral when it is combined with 
the 5-terms appearing in the definition of D; however, these contributions cancel. The 
remaining quantities depending on L^> are found in D\ and are already O(At), so 



D 



1 TAt K" , 
+ —At 2 

2 2 16 
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they provide contributions of 0(At 2 ) when integrated over t w < t' < t, if the remaining 
quantities in the integral are 0(1). Again, this can only be realized in the single integral 
by setting the three time function CC to its zero order value, and in the coupling with 
the 5-terms in the double integral, but these terms cancel. Also, in the final result there 
is a leading order cancellation of the terms proportional to K"(0) and one obtains 

V c (At) = 2T 2 At 2 (4.30) 

This agrees in both the At-dependence and the prefactor with the result of the high-T 
expansion. It is worth noting that in order to get the short time behaviour of the 
correlation variance, we have expanded CCit' — t' w ) around CC(0). This quantity 
diverges in d < 4 at T = T c . But the fact that it cancels from the leading short- 
time behaviour should mean that (14.301) remains valid: in principle we just need to 
regularize in some way, e.g. by keeping t w large but finite, then perform the short-time 
expansion and finally remove the regularization, which should lead back to (I4.30p . 
For the susceptibility variance, as given in (13.431) . we can expand D x as 

D x = At 5{t - t') - 6(t' - t w ) + <D(At)6(t' - t w ) + <D(At 2 )6(t w - t') (4.31) 

The first two terms and the remainder give contributions of O(At) and 0(At 2 ) 
respectively to the integral over t' . Keeping only the first two terms and replacing 
CC{t' - tj by CC(0) in (E3D should then give the leading term in V x of £>(At 2 ); but 
this cancels because J dt'[At5(t — t') — 6(t' — £ w )] = 0. The same argument shows that 
to 0(At 3 ) the cross terms between the O(At) and 0(At 2 ) contributions from ( 14.311) 
cancel. The only remaining 0(At 3 ) term is then 

V x = \f dt'dt'JAtS(t - t') - 9(t' - tJ][At6(t - t'J - 9(t' w - t w ))CC(t' - t'J (4.32) 

= ^At 2 CC(0) - At f dt' CC(t - t') + f dt' f dt' w CC(t' - t'J (4.33) 

One now expands CC(t' - t'J = CC(0) - 2T(f - t'J for t! > tj using (gSD, to find 
that the At 2 CC(0)-terms cancel as expected, leaving 

V x (At) = -T 3 At 3 (4.34) 
3 

This again agrees with the high-T expansion in both the scaling with At and the 
prefactor. 

For the covariance, finally, performing the perturbation expansion in small At shows 
that the leading short-time contribution is only of fourth order in At. The temperature 
dependence of the prefactor is quite complicated and involves the kernel L^ 2 \ We only 
show here the limit of the prefactor for large T, which is 

AdT 2 

V Cx {At) = -^-At 4 (4.35) 

in agreement with (14.221) . Correspondingly, the correlation coefficient is again negative, 
growing in modulus initially as 7 ~ —At 3 / 2 . The fluctuation slope Xr can also be 
expanded as in the high-T case ( 14.241) . leading again to a quadratic short-time increase 

A fl ~ ~ At 2 (4.36) 

Vc 
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4-4- Large time differences 

We next turn to the behaviour of the correlation and susceptibility fluctuations at 
large time differences, at equilibrium at temperatures above criticality; the dynamics 
at criticality exhibits qualitative differences and is considered separately in the next 
section. The asymptotic behaviour of (the equilibrium forms of) K, and CC 
is an exponential decay. For K and CC this follows directly from (14.41) and (14. 8p . 
where all Fourier modes decay as exp(— 2z cq At) or faster. The Laplace transform 
K{s) then has all its singularities bounded away from s = 0, and the same follows 
for L^{s) from (14111 . Thus, looking at D(t,t w ,t') as given by (fl~T3l) . the K{At/2)- 
prefactor ensures that the 5- contributions to D decay exponentially for large time 
differences, and the same is true for the Dx-term found in (14.141) . The D 2 -term 
given in ( 14.151) has the same behaviour. This is obvious for the first term; for the 
second term, bounding both K[t) and Lfi>{t) by exp(— ct) shows that the integral is 
bounded by Atexp[— c(t — t')] < At exp(— cAt/2). Thus, for large At, all the non- 
Gaussian corrections in (13.411) . as well as the two-time Gaussian term CC(t,t w ), decay 
exponentially to zero; the asymptotic value of the correlation variance is then given 
by the time- independent Gaussian term J(dq) Cq(t,t)C n (t w ,t w ) = J(dq) T 2 /(u + z cq ) 2 . 
This increases as the temperature is reduced towards T c ; the limit value for T — > T c 
is finite for d > 4 but infinite for d < 4. From the reasoning above it follows that for 
T > T c the approach of the correlation variance to its asymptotic value for At — * oo is 
exponential in At, up to power law factors. 

Analogous reasoning for D x in ( 14.161) leads one to discard as subleading for large 
At the D x term in (14.181) and the terms proportional to K(At/2) appearing in (14.171) 
for D\ . The asymptotic susceptibility variance can therefore be found from (13.431) by 
replacing D x with 

QL* = S(t - + \K'{*f) + 2TK(tf) - jT dt" LW(t» - 0)K(t=f) (4.37) 

We have added the subscript "short" because all contributions retained in -D^ hort decay 
exponentially as t — t' increases, concentrating the "mass" of the integrals over t' and t' w 
in (gaUD into the regions t-t' = 0(1) and t-t' w = 0(1). (The factor CCit! - t' w ) does 
not affect this reasoning as its values are also largest when t' is close to t' w .) Nevertheless, 
because CCit 1 — t' w ) does vary significantly on 0(1) timescales one cannot simplify the 
expression for the asymptotic value of V x further, beyond the replacement of D x by 
D* hort in (13.431) . Barring accidental cancellations (which, by continuity with the nonzero 
result for T — ► oo, one does not expect), the result will be nonzero. The approach to 
the limit will again be exponential in At. 

Finally, in the cross-correlation (13.441) there are no 0(1) contributions that survive 
in the long-time limit, because every term is proportional to D and thus decays 
exponentially with At. Qualitatively, then, the behaviour at finite T > T c is the same 
as that for T — ► oo, with the variances of correlation and susceptibility approaching 
nonzero asymptotic values exponentially fast in At, and the covariance decaying to zero 
(from below) in the same manner. 
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5. Quenches to T c , d > 4 

Above we derived the large- N statistics of the correlation and susceptibility fluctuations 
for quenches to temperatures above criticality. In this section and in Sec. [6] we consider 
quenches directly to criticality, so from now on T = T c = (J(dq) 1/uj)" 1 . In principle we 
then expect aging effects [2J [33], but it will turn out that for the fluctuation statistics 
these are largely negligible as long as we are in dimension d > 4. We therefore consider 
first the situation in equilibrium at criticality in d > 4, focussing on large time differences 
At; the short-time limit does not need to be analysed again here because the results in 
Sec. 14.31 apply even at T = T c . As discussed in more detail in Sec. El for d < 4 one has 
to keep t w finite to avoid the appearance of infinite terms, i.e. one has to look directly 
at the non-equilibrium situation. 



5.1. Equilibrium 

The asymptotic behaviour of the equilibrium forms of K, CC and for T — > T c is 
quite different from the high temperature phase b vanishes. In the equilibrium 

form (14.41) of the kernel K, the integral is for large time-differences dominated by small 
uj. Because uj « q 2 for small q = |q|, the phase space factor in the q-integrals is 
(dq) = adduj uj d l 2 ~ 1 for small q or uj, with the proportionality constant 

a d ={^)- d l 2 V-\d/2) (5.1) 

Then from ( 13.151) one finds for large At 

K(At) =a d [ dujuj^ 2 -e- 2ujAt = k d At^' 2 (5.2) 

J UJ 

and correspondingly, from the small- s expansion of (14. lip . 
rffl/A , x f At^-V/ 2 ford>4 

L»(*) = A<{ M(d _ m fmJ<4 (5.3) 

with k d and ^-dependent constants [33] . 

In equilibrium at criticality the function CC(At) also decays as a power law for 
At ^> 1. This can easily be worked out in d > 4: 

(dq) —e~ 2ulAt = 2T / d£ K(t') ~ At (4 - rf)/2 (5.4) 

The three-time function CC(t,t w ,t') given in (I4.10p increases with £' up to t' = t w ; in 
the range t' — t w . . . t it then remains constant and equal to CC(At/2). 

For the correlation variance (13.411) . the first of the two Gaussian terms is a constant 
of order unity, J(dq) C q (t, i)C q (i w , t w ) = j(dq)T 2 /u 2 . The second one, CC(At), is 
proportional to At^ 4 "^/ 2 for large At from (15.4H . One can show that these are in fact 
the two leading terms, with the non-Gaussian corrections contributing at most 0(At 4 ~ d ) 
asymptotically. (See [42] for a detailed discussion of the scaling of these terms.) Note 
that because the prefactor of CC(At) ~ At^ 4- ^/ 2 in (13.411) is positive, Vc approaches 
its (positive) asymptotic value from above; given that it starts at zero at At = 0, 
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it is therefore non-monotonic in At. Because Vc-(At) must depend continuously on 
temperature, at least at finite At, this non-monotonicity must then be present also in a 
range of temperatures above T c . 

To understand the asymptotics of the susceptibility variance (13.431) it is useful to 
decompose 

D*(t, t w , = Dl OTt (t, t w , t') + Dl ng (t, t w , t') (5.5) 

where 

Dl ng (t, t w , f ) = -K(f) 8(t - W*> *w, t')6(t' - t w )+D*(t, t w , t')9(t w - f) (5.6) 
and 

£>W*> t w , t') = -K(f) [2T - J* alt" - f )] (5.7) 

while D x is as written in ( 14. 18ft . The short-time part of (15.51) . defined in (14.37p . decays 
on timescales t — t' = 0(1) and its integral over t' is of order unity. More precisely, 
Laplace transforming D* hoTt (t , t w , t') w.r.t. t — t' and expanding for small s gives 

Dl oA s) = 2K(2s)[s + 2T - L W ] = « 2 ^p^S^ ^ 

~ 2 + 2-^ S ( rf - 4 )/ 2 (l - 2( d - 4 )/ 2 ) (5.9) 
K(0) v ; V 7 

where we used (14.111) and the small s-expansion of K(s) = K(0) — cs^ -4 -*/ 2 [33]. This 
shows that the integral of D* hort (t — t') over all t — t! equals 2, and that D* hort (t — t') 
decays as (t — t')( 2 ~ d )/ 2 for large At. The complementary long-time part D* ong , on the 
other hand, can be seen to have a structure similar to that of D: the 5-term has weight 
~ At (2 ~ d)/2 , -Df jlong is of order At (2_a!)/2 (and constant for t — t' > 1), and D\ is of order 
L^(t — t') ~ (t — t'Y 2 ~ d ^ 2 . Decomposing now the product D X D X in (13.431) according 
to (I5.5p . the -D^hort-^short t erm gives the leading asymptotic term of V x , a constant of 
order unity. We omit here the detailed analysis of the scaling of the other terms, which 
can be found in [22] , and point out only that overall V x approaches its asymptotic value 
via leading power law terms ~ At^ 2- ^/ 2 and ~ At 4_d ; the former dominates for d > 6, 
the latter for 4 < d < 6. There appears to be no simple way of estimating the sign of 
these terms to verify whether V x , like Vc, has a non-monotonic At-dependence. 

In the co variance (13. 44ft there are, as in the case of T > T c , no 0(1) terms that 
survive for long times. To obtain the leading decay to zero, one considers first the single 
integral involving D x and the three-time function CC(t, t w , t'). The leading contribution 
arises from the short-time part D* hort of D x , giving w —2CC(t,t w ,t) = — 2CC{t — t w ) ~ 
— At*- 4-6 '-'/ 2 , so that Vcx(At) approaches from below as — At^ A ~ d ^ 2 . The scaling of the 
remaining subleading terms is analysed in [42] . 

A schematic plot of the full At-dependence of the variances and the covariance at 
criticality and d > 4 is shown in Fig. [51 The non-monotonicity of the correlation variance 
(and, possibly, the susceptibility variance) is consistent with a scenario in which for both 
short and long time regimes local correlations and responses display minor fluctuations 
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Figure 5. Schematic plot of the Ai-dependence of the variances, Vc and V x , and 
covariance, Vc x , for critical coarsening and d > 4. These functions increase from to 
0(1) values as power laws of At given respectively in (|4.30|1 . (|4.34p and (|4.35|) . and then 
decay asymptotically with negative powers of At. Vc and Vc x are non-monotonic; for 
V x it is not obvious whether the asymptotic value is approached from above or below. 

across the system and dynamical histories, whereas for intermediate time they depend 
strongly on the particular state reached by the system in a given time interval, making 
fluctuations across dynamical histories large. This type of behaviour is quite generic: 
indeed, as explained after (|2.4p . Vc is the four-point function commonly used to quantify 
dynamic heterogeneities, and this is found to exhibit a maximum as a function of At 
for many systems with slow dynamics [TOj HH [H]. Note though that typically (e.g. in 
coarsening below T c [231 [21]) the position of the maximum scales with the age of the 
system whereas here, for coarsening at criticality above the upper critical dimension, it 
occurs for an age-independent time difference At. 

Gathering the above results for At -> oo, V c (At) -> 0(1), V x (At) -> 0(1) and 
Vc x (At) — > —At( 4 ~ d ^ 2 , yields for the correlation coefficient the power law decrease 
j ~ — At^ d ^ 2 . The FDR for the local fluctuations Xq as given in (13. 5 7ft is then positive, 
but to say more one would need to know whether Vc and V x is bigger for At — > oo. If 
Vc > V x , the elliptical contour of P(C, T\) has its main axis along the C direction and, 
as one can show formally by expanding (13.571) for small 7, X^ ~ —7 ~ 

At (4-d)/2 decays 

to zero. Conversely, if V x > Vc the main axis of the ellipse is along the x direction 
and the fluctuation slope Xr ~ — 1/7 ~ At^ d ~ A ^ 2 becomes vertical for At — >• 00. In 
the limiting case where Vc and V x are equal asymptotically, (V^/KJ 1 / 2 — (V x /Vc) 1 ^ 2 
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in (13.571) would decay as At ( - 4 ~ d ^ 2 , being controlled by the leading corrections to Vc] 
because this is proportional to 7, a finite asymptotic value of would result. From 
the high-T expansion of Sec. 14.21 one can see that at order 1/T 2 , differences between the 
asymptotic values of Vc and V x appear. These suggest that the scenario that should 
apply is the one where V x > Vc asymptotically and therefore the main axis of the 
fluctuation ellipse is along the x direction. 



5.2. N on- equilibrium, d > 4 



Now we study the behaviour of correlation and susceptibility fluctuations for the genuine 
out-of-equilibrium dynamics after a quench to criticality, focussing on d > 4 as in the 
previous subsection. Here we get correction factors with respect to the equilibrium case 
which become important in the aging regime, t — t w ~ t w . More precisely, aging effects 
appear in the long time (t, t w ^> 1) behaviour of two-time functions via scaling functions 
of the time ratio t/t w that modulate the equilibrium part, e.g. 

(5.10) 
(5.11) 



K(t,t w ) = K cq (t-t w )F K 
L( 2 )(t,t w )=Lg)(t-t w )^ L (f 



with Fk(X) = Fl(X) = 1 [33J. Here and throughout the remainder of the paper we 
distinguish the TTI equilibrium contributions with the subscript "eq". The two-time 
correlation C(t,t w ) has a similar scaling behaviour as can be seen by expressing it in 
terms of the kernel K: 



C(t,t w ) 



(dq) -Rq(t, t w )Cq(^w, ^w) 

g(t) 



J^J{dq)e~^C (l {t W) t w ) 



-K(i,t w ) =K cq (f ) 



git) 



Mi) 



(5.12) 
(5.13) 



'g{t)g{U 

The long-time behaviour of the function D from (13.401) that defines the correlation 
variance is thus given by (using (13.471) . (I3.48P and (13.501) ) 

g(t) 



D(t,t w ,t') = |*: eq (f ) 9{t) Fk(*¥) [6(t - t') + 5{U - 0] 



'g(t)g(t v 

+ D^t, t w , t')6(t' - t w ) + D 2 (t, t w , t')6(t w - t') 
where x = t/t w and 



D 1 (t,t w ,t') = -K eq (f) 



g(t) 



and 

D 2 (t,t w ,t') 



-K ^ 

2 eq \ 2 



g{t)g{t w ) 
g(t) 



K 



x+1 
2 



2T-J t t i drL^(T-t')T L ^ 



(5.14) 



(5.15) 



g{t)g{U 



1 N ' Jt^ x 



+ 



git) 



g(t)g(t v 



f drK eq (i- r) ^(i) Lg(r - t)F L {$ 

J tvj 



(5.16) 
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For d > 4 and long times g(t) ~ const and all the factors involving g(t) in the above 
expressions can be dropped. 

In order to get the long-time expression for D x from (13.521) we need the non- 
equilibrium form of x(Mw)- In equilibrium, 1 — Tx eq (A£) = K cq (At/2) and so we 
consider the same combination here: 



1 - T X (Mw) = 1 - T I (dq) f dt' e'^'K 

J Jt w \ 



g(t) 



(5.17) 



This can be rewritten, by adding and subtracting the quantity T J(dq) // dt' e ^ 1 ^ 
1 - K cq (At/2), as 



1 - T X (t, t w ) = K cq (f ) - T f (dq) f dt' <?W 



M 9(t) 



(5.18) 



and extracting a factor of K eq (At/2) from the second term yields 

l-T X (t,t w )=K cq (f)F y (x" 
Bearing in mind that / (dq) e^^i^ 1 ') - 



(5.19) 

K' (t — t')/(2T) the aging function T x can be 



written as 



TJx) = 1 + 



9(f) 

<?(*) 



(5.20) 



2 K eq (At/2) 

Because g(t) approaches a constant for large times in d > 4, the integral term in fact 
vanishes in the limit and there is no aging correction: T x (x) = 1. Inserting (15.191) 
into (I3.52p . (13.531) and (13.541) one finds the long-time non-equilibrium form of D x 



D*(t,t w ,t/) 
with 

D x (t,t w ,t') 



l-K cq (f)F x (x)} 5(t-t') 
+D x (t, t w , t')9(t' - t w ) +D x (t, t w , t')9(t w - t') 

+ 2T [K^(*£) Fxi?) ~ K«i(¥) *x 
- tdfL^f'-t')^) [Kj^)r x {$)-Kjtf)F x 



(5.21) 



(5.22) 



and 

D x (t, t w , t 



t-t' 



• r x\t" 



K eq (f)f x (x)] (5.23) 
(5.24) 



dfL$( t »-t>)F L Ci] 

It will be useful to separate short and long time parts in D x again: 

D*{t, t w , = D x hoit (t, t w , if) + D x ng (t, t w , t') 

We arrange the terms so that the first term is identical to its equilibrium 
counterpart (14.371) except for a slowly varying aging correction, giving 

D x lort (t,t w ,t') = F x (±) \s(t - + \K q (±f) + 2TK cq (^f 
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dt"L^{t" - t>)K4^)]p{i! - t w ) (5.25) 

D? ong (t,t w ,t>) = -K(f)f x (x)5(t-t') + Dl long (t,t w ,t>)6(t f -t w ) 

+DZ(t,t w ,t')0(t w -t') (5.26) 

^W*> W) = ^ eq (^)^(|) ~K^(f )T x (x)[2T-£dt"L^(t"-l/)T L (f) 

~ f t ,dt"L%{t" - t>)K cq (^f) [t l (^) T x ([±) - T x ^)] (5.27) 

We have kept the factors of T x so that the expressions are valid also in d < 4, for later 
use. For our current case (d > 4), one has JF X = 1 and JF^ = 0. 

To deduce the behaviour of Vc, V x and Vc x we finally need the aging corrections 
to the function CC defined in (13.371) . (I3.38[) . One uses (I3.10P and the long-time scaling 
of the equal-time correlator C q (t w ,t w ) = (T /co)Tc(ujt w ) with Tc{w) — > 1 for w — > oo 
and J-c( w ) ~ 2iy for w — > 0; in d > 4, JFc(ty) = 1 — exp(— 2w) [33]. The ratio of the 
three-time function to its equilibrium counterpart is then, for t' > t w , 

CC(t,t w ,t>) J(dq) (T 2 /m 2 )F c (ut')F c (Lut w )e-^-^ 
CC eq (t,t w ,t>) f(dq){T*/u2)e-»(*-*~) { - ) 

In the aging regime where t — t w ^> 1 and so only w € 1 contributes one can replace 
(dq) ~ did uj d l 2 ~ 1 ] rescaling lo to w = ujt w then gives a function of x = t/t w and y = t' /t w : 

CC(t,t w ,t') J dw w( d -V/ 2 F c (wy)Fc(w)e-^^ 

7- = Fcc(x,y), T cc (x,y) = F^T^ev^fa^ ( 5 - 29 ) 



CC eq (t,t w ,f) oov oov ' yy J^^( rf -6)/2 e --(-i) 

With the same approach one finds for t' < t w , i.e. y < 1, 

_ Jrf WW ( rf - 6 )/ 2 ^( W i/)e-^ +1 -^ 

The function Tec a ls° governs the aging corrections for the two-time CC, according to 

CC(t, t w ) CC(t, t, t w ) <p cc ^i j/^j ^ 3^ 

CC eq (t — t w ) CC eq (t, t, t w ) 

We can now proceed to study what, if any, aging corrections there are for the 
correlation and susceptibility (co-)variances in d > 4. The scaling functions Tk, Tl, 
JF X , Tec are an bounded by 1. One can check that the presence of the aging corrections 
in D and D x does not increase the order of the various terms and that the asymptotic 
values of Vc and V x are identical to the ones in equilibrium, with aging corrections 
visible only very weakly in the decay to these asymptotes. A detailed discussion of 
the effects of aging corrections on Vc and V x is provided in [42]. Here we mention 
only that the leading correction to the asymptotic value of Vc displays a crossover 
from the value CC eq (t — t w ), obtained in the near equilibrium regime x ~ 1 (where 
Tcc(^-A/ X ) ~ I); to an x-independent negative constant of order t ( £~ d ^ 2 for x 3> 1. 
This negative contribution provides a small downward shift in the value of Vc that is 
reached asymptotically, as t becomes large (t ^> t w ) at fixed t w . 

The only significant aging effect in critical coarsening for d > 4 is 
visible in the covariance, whose leading term can be shown to be [42] Vc x = 
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-2CC cq (At/2)F cc (x,x) 4 4 " d)/2 (x - iy 4 ~W 2 F C c( x , x ), that is the equilibrium 

result supplemented with an aging correction Fccfa, x )- F° r x ~ 1 the latter equals 
one as it should to reproduce the equilibrium result; for large x, one finds from (15.291) 
that J r cc( x ,y) ~ vl x<1 an d so J r cc( x , x ) ~ \ jx. 

6. Quenches to T c , d < 4 

In this section we will study the fluctuations in the out-of-equilibrium dynamics of 
the spherical ferromagnet quenched to criticality in dimension d < 4. Here we 
cannot start from an equilibrium calculation and later account for aging corrections 
because a naive equilibrium limit leads to the appearance of infinite terms: in the 
correlation variance (I3.4ip . for example, the first Gaussian term becomes J(dq)T 2 /uj 2 
for t, t w — > oo which is infinite in d < 4. Thus, we need to look directly at the non- 
equilibrium situation. Specifically, we will consider the aging limit t, t w — > oo but at 
fixed x = t/t w > 1. We omit the calculations for this scenario as they are somewhat 
technical; the interested reader can find full derivations in Ref. [42J. 

Our analysis shows that in the aging regime the variances and the covariance all 
scale as t^~ d ^ 2 times a function of x. The dependence on x implies that the relevant 
timescale on which the (co-)variances vary is At ~ t w , in contrast to the case d > 4 where 
this timescale is At = 0(1) independently of the age t w . In this sense the behaviour 
for d < 4 is similar to what is seen for coarsening (in general d) at T < T c [23] [24]. 
Interestingly, however, the amplitude of the (co-)variances grows with the age t w only 
as t^~ d ^ 2 at criticality, whereas below T c it scales (at least for Vc) in the naive way as 
the domain volume ~ (t\[ 2 ) d . 

The dependence on x of the (co-)variances can be evaluated numerically from the 
scaling expressions for the aging limit; the results are shown in Fig. [6] for d = 3. We 
also study numerically the x-dependence of the resulting correlation coefficient 7 and 
the fluctuation slope Xg_, as displayed in Fig. [7J In both of these quantities the t^~ d ^ 2 
prefactor from the (co-)variances cancels, so they depend solely on x. 

One can analyse the behaviour of the quantities above in more detail for the two 
opposite extremes e = 2; - 1 C 1 and x ^> 1. In the former case we expect to recover 
quasi-equilibrium behaviour with all dependences being only on time differences. One 
can show by an expansion in e [42] that indeed the (co-)variances all grow as e^ d ^ 2 
to leading order; combining this with the overall t^~ d ^ 2 scaling, one gets a TTI time 
dependence as expected, proportional to (t — £ w )( 4_d " 2 . The prefactors Uc, v x and vq x 
of this power law increase in Vc, V x and Vc x are plotted as functions of dimensionality 
in Fig. El Given that Vc, V x , Vc x all have the same scaling with t — t w in this regime, the 
correlation coefficient and fluctuation slope also have nontrivial values; these are plotted 
against d in Fig. [91 Note that because we are approaching the TTI regime via the limit 
x — > 1 of an aging calculation, where t — t w ~ t w ^> 1 always, these results are valid for 
t — t w ^> 1 only. For t — t w <C 1, on the other hand, the results of Sec. 14.31 will apply and 
both 7 and X& will tend to zero in the limit t — > t w . If we were to plot the contour lines 



Dynamic heterogeneities in critical coarsening 33 




0.01 0.1 1 10 100 1000 10000 
x-1 



Figure 6. Log- log plots of correlation and susceptibility variances Vc (top left), 
V x (top right) and covariance Vc x (bottom) versus x — 1 = (t — t w )/t w , for critical 
coarsening in dimension d = 3. As in all following figures, the (co-)variances have been 
divided by tw d ^ 2 to get functions of x only. The theoretically expected power laws 
for large x are given, respectively, by Vc ~ x°, V x ~ x^ d '' 2 and Vc x ~ — x 2 - 3d / 4 
and are represented by the dotted lines on the right of each graph. The 0(1) value 
approached by the correlation variance is represented by the Gaussian term and it 
can be calculated analytically this value is represented by the horizontal dotted 
line in the first graph. The initial increase ~ (x — l)( 4 ~ d )/ 2 for all three quantities is 
indicated similarly on the left. Note that for V x the initial and asymptotic power laws 
are identical, so that the log-log plot overall is close to linear. 

of the distribution of (C, T\) on an FD plot the initial section (bottom right) would 
therefore look similar to Fig. HJ but then the ellipses would grow as (t — t w )^ 4_d ^ 2 as the 
top left hand corner of the plot is approached and their principal axis would approach 
a limiting slope given by in Fig. [9] as a function of d. The genuine aging effects 
occurring for x > 1 would not be visible because they are all compressed into the top 
left corner of the plot in the limit t w — ► oo. 

Analytically one can obtain relatively simple expressions for the prefactors uc, v x , 
vc x in the limits d — > 2 and d — > 4 [12]. We mention here only that they all vanish as 
power laws of d — 2 in the limit d —*■ 2. This makes sense intuitively in that for d < 2 no 
phase ordering takes place and so the fluctuations arising from the coarsening dynamics 
should vanish as d — >• 2 from above. 
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Figure 7. Log-log plot of the absolute value of the correlation coefficient 7 (left) and 
the fluctuation slope Xq (right) versus x — 1 = (t — t w )/t w for critical coarsening in 
d = 3. The modulus of the correlation coefficient 7 is always smaller than 1, as it 
should be. For large x, 7 decays as — x( 2 ~ d ^ 2 , as shown by the dotted line on the right 
of the plot. The separately calculated limit value for x — > 1 is indicated by the arrow 
on the j/-axis and is certainly plausible as an asymptote for In (x — 1) — > —00 of our 
numerics for x close to 1. The fluctuation slope Xfl behaves for large d ' i . This 

predicted power law is represented by the dotted line on the right, while the value that 
should be approached for x — + 1 is indicated by the arrow on the y-axis. 





Figure 8. Left: The prefactors of the (co-)variances in the TTI regime x — t/t w w 1, 
where they are proportional to (< — t w )^ 4_d '/ 2 , are plotted versus d. All three prefactors 
diverge as 1/(4 — d) for d close to 4 and vanish as power laws for d near 2, specifically 
vc ~ id— 2) 2 and v x ~ ~ (^~2) 3 . Right: Plot of the prefactors normalised by their 
predicted limiting behaviour as 5 = (4-d)/2 42J, 167r 2 <5 i / c /(3T c 2 ), 4tt 2 5v x /T* and 
—8ir 2 5vc x /Tc, versus d The inset shows the correlation variance prefactor normalized 
by its predicted limiting behaviour as 5' = (d — 2) /2 — > 0. 
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d d 

Figure 9. Correlation coefficient 7 and fluctuation slope Xa for x = t/t w — > 1, plotted 
versus d. The circles show the limits of 7 and X^ for d — > 4, calculated analytically and 
given respectively by — 1/ V3 and (l + VT7)/4 (see 42 ) . The inset on the left shows the 
correlation coefficient normalized by its predicted power law as it approaches d = 2. 
The loss of accuracy around d = 2 is due to numerical inaccuracies in the evaluation 
of the correlation variance as d gets very close to 2 [42] . 



Turning to the opposite limit of large x, we find [42] that the correlation variance 
Vc is dominated by the Gaussian term asymptotically, matching quantitatively the 
behaviour in the regime d > 4. Here, however, the first subleading correction to the 
constant asymptote is negative, so that the approach is from below, in contrast to d > 4 
(see Fig. [6]). Quantitatively the correction term is small already in d = 3 and numerical 
evaluation shows that is gets progressively smaller as d increases to 4. For the response 
variance and the covariance we get V x ~ x( 4_c 0/ 2 and Vc x ~ —x 2 ~ zd l A . All of these 
scalings, as well as those for x — > 1, are in agreement with our numerical evaluations as 
shown in Fig. [6] above. 

We discuss briefly the consequences of the above results for the large-x behaviour 
of the contour ellipses of the joint distribution P(C, T%) of the fluctuating correlation 
and susceptibility. Firstly, due to the different scaling with x of the correlation and 
susceptibility variances, as x° and x^~ d ^ 2 respectively, these ellipes become increasingly 
elongated in the susceptibility direction as x grows large. Using also the scaling of the 
covariance Vc x ~ — x 2 ~ M//4 , the correlation coefficient 7 from (13.551) decays to zero 
aymptotically as — x^ 2 ~ d ^ 2 . Finally we consider the fluctuation slope Xr. Looking 
at (13371) . one can check [42\ that this is given by X fl = -2(V x /V c ) 1/2 /{2-f) = -V x /V Cx ~ 
x d ^. The large-x divergence of Xr as x d//4 for any 2 < d < 4 is consistent with the fact 
that the joint distribution of correlation and response fluctuations grows more quickly 
in the susceptibility direction than along the correlation axis. It also matches our finite 
x numerics, see Fig. [7J 
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7. Discussion 

We have analysed fluctuations in the coarsening dynamics of the spherical ferromagnet 
after a quench, specifically the leading 1/yN fluctuations of local correlations and 
susceptibilities spatially coarse-grained across the entire system. Our work was inspired 
by general theories regarding the nature of correlation and susceptibility fluctuations 
in aging systems [HJ [T7J [TH [19] . Our study significantly extends the scope of previous 
(zero-temperature) calculations of correlation fluctuations in the spherical model [25] by 
keeping track of non-Gaussian fluctuations. This enables us to calculate explicitly the 
susceptibility fluctuations, which in a Gaussian approximation would vanish identically. 
The nature of our approach, which treats the non-Gaussian effects perturbatively, means 
that we cannot analyse quenches to below the critical temperature; however, we can 
access the interesting regime of coarsening at criticality, where our results are the first 
of their kind. 

We discussed carefully in Sec. |2] possible definitions of coarse-grained fluctuating 
correlations C and responses \. It turns out that for the fluctuation statistics (in non- 
disordered systems such as the one studied here) it does matter whether the underlying 
local functions are measured directly, or indirectly via quenched amplitudes that define a 
randomly staggered magnetization observable: only the former choice gives correlation 
and susceptibility variances that scale in the same way with system size N. These 
considerations should be of general relevance to other systems where a coarse-graining of 
the local correlation and response across an entire finite-sized system is desired. Coarse- 
graining over smaller volumes does not produce interesting results in the spherical model 
because the susceptibility fluctuations are small (~ iV~ 1//2 ) but correlated across the 
entire system. 

In Sec. [3] we used the 1/vJV expansion of the non-Gaussian fluctuations [33] to 
derive general expressions for the correlation and susceptibility variances and covariance. 
These are exact to leading order in 1/N, where the joint distribution of C and x is 
Gaussian. In addition to Vc, V x and the covariance Vc x , this distribution can be 
characterized by the correlation coefficient 7 (see equation (13.551) ) and the negative slope 
of the principal axis of the elliptical equi-probability contours, Xr (see equation (13.57P ). 
The definition of Xq was chosen such that, if the predictions for glassy systems (such 
as spin glasses) with a global time reparameterization invariance [161 (HI EEl EE] applied 
also to coarsening systems, X& should be close to the fluctuation-dissipation ratio (FDR) 
X that relates the variations with time of the average susceptibility and correlation. 

In Sec. H]we considered first quenches to T > T c , where after fast initial transients 
the dynamics is in equilibrium. Analytical results were obtained in the limit of high 
temperatures: here the probability contours rotate with increasing time difference 
At = t — t w from a horizontal orientation (Xr = 0) to a vertical one (Xr = 00). 
At the same time the correlations between correlation and susceptibility fluctuations 
become weaker and weaker: the contours become approximately circular, showing an 
effect opposite to the progressive narrowing of the contours around the slope of the 
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fluctuation dissipation plot (X fl = 1 for large T) that would be expected for spin glasses 
and similar systems [161 El HI EES] • Qualitatively this behaviour remains the same also 
for quenches to lower temperatures above T c ; the results for small time differences At, in 
particular, depend on T only through prefactors. The correlation coefficient 7 between 
the fluctuations of correlations and susceptibilities is always negative, corresponding to 
a positive fluctuation slope X fl . 

The more interesting quenches to criticality were studied in Sec. [5] (for dimension 
d > 4) and Sec. [6] (for d < 4). In the former case, we found that out-of-equilibrium 
effects are weak and one can directly analyse the equilibrium dynamics. Interestingly, 
the correlation variance - which is identical to the four-point correlation function often 
used to characterize dynamic heterogeneities - displays a maximum as a function of At, 
suggesting as in other glassy systems that there is a well-defined timescale on which 
fluctuations between different dynamical trajectories of a system are largest. However, 
even though the coarsening dynamics has a growing lengthscale that increases with the 
system age t w in the standard diffusive manner, £(t w ) ~ tl/ 2 , neither the timescale of 
the maximum in Vc nor its amplitude change with age. This is in contrast to the case 
of coarsening below T c , where timescales grow with the age and the variance has the 
natural scaling with £ d (t w ) ~ t^ 2 |23j. 

Below d < 4 one has to look directly at the non-equilibrium situation: a naive 
equilibrium limit yields infinities that need to be regularized by initially keeping the 
age t w finite. We found that Vc, V x and Vc x a ^ sca le as tw~ d ^ 2 times functions of the 
time ratio x = t/t w . Looking at the details of the x-dependence, we saw that in the 
quasi-equilibrium regime x ~ 1 time-translation invariance is restored as expected, with 
all (co-)variances scaling as t^~ d ^ 2 (x - l)( 4 - d )/ 2 = (t - t w )( 4 "^/ 2 for t - t w » 1. Unlike 
the case d > 4, the correlation coefficient is finite in this regime, but the corresponding 
fluctuation slope does not seem to be related to the FDR - which is X = 1 at 
quasi-equilibrium - in any simple way. In particular, Xq grows monotonically from a 
vanishing value at d = 2 to the non-trivial limit (1 + -\/l7)/4 in d = 4. In the genuine 
aging behaviour that follows for larger x the correlation coefficient decays to zero and the 
fluctuation slope X fl diverges towards large positive values as x grows, both reflecting 
the progressive stretching of the probability contours along the susceptibility axis. 

From a more general point of view, our results show clearly that heterogeneities 
are present in coarsening at criticality above the upper critical dimension, as detected 
e.g. via maxima in Vc in d > 4. However, the lack of a dependence on age t w in the 
relevant timescales (At = 0(1)) and amplitudes (Vc = 0(1)) is somewhat surprising. 
Interestingly, the maximum in Vc is seen to disappear below the critical dimension. The 
timescale on which Vc varies then has a conventional aging form (At = 0(t w )) while 
its amplitude Vc = 0(t^~ d ^ 2 ) is not related in any obvious manner to the growing 
correlation volume of order t^ 2 . 

It will be an interesting challenge to see whether the general features of our results, 
and in particular the different t w -scalings of the correlation variance above and below 
d = 4, can be understood from general scaling or field theoretical approaches to critical 
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coarsening [28J. One would also like to extend our considerations to genuinely short- 
ranged systems: the spherical model is somewhat unusual in that the spherical constraint 
generates a weak but long-range interaction. The 0(n) model in the limit of large n may 
be a suitable candidate here; preliminary work suggests that much of our perturbative 
approach for analysing non-Gaussian fluctuation effects would transfer to this scenario. 
Finally, it is clear from our results that fluctuations in critical coarsening display very 
rich behaviour that cannot simply be deduced from the properties of the average 
fluctuation-dissipation relations, and it remains to be seen whether alternative ways 
can be found of rationalizing the kind of effects thrown up by our exact calculations. 
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